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ABSTRACT 

The Laser Interferometer Space Antenna mission is a planned gravitational wave detector consisting of three spacecraft 
in heliocentric orbit. Laser interferometry is used to measure distance fluctuations between test masses aboard each 
spacecraft to the picometer level over a 5 million kilometer separation. Laser frequency fluctuations must be suppressed 
in order to meet the measurement requirements. Arm-locking, a technique that uses the constellation of spacecraft as a 
frequency reference, is a proposed method for stabilizing the laser frequency. We consider the problem of arm-locking 
using classical optimal control theory and find that our designs satisfy the LISA requirements. 

INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) is a space-based gravitational wave detector under joint development by 
NASA and ESA [1-3]. The goal is to measure gravitational waves, ripples in the curvature of space-time produced by 
astrophysical events such as black hole collisions. Gravitational waves produce an oscillating tidal strain in space, the 
effect of which is to modulate the distance between a set of separated inertial masses. The amplitude of the oscillations is 
proportional to the nominal distance between the masses. Therefore the wave’s amplitudes are characterized by the 
strain, or ratio of change in length to nominal length. For astrophysical sources of interest to LISA, typical strain 
amplitudes are 1 0' 21 and have Fourier frequencies between 10' 4 Hz and 10' 1 Hz. 

The LISA instrument consists of atriangular constellation of three spacecraft in heliocentric orbit at 1 AU, trailing 
approximately 20 degrees behind Earth. Each arm in the constellation is nominally 5 million km long and varies by as 
much as one percent with an annual period. Laser interferometry is used to monitor distance fluctuations between freely 
floating masses that are contained within the spacecraft at each end of the arm and are isolated from external 
disturbances using the technique of drag-free control [4,5]. A measurement resolution of approximately 10 pm (1 pm = 

1 0' 12 m) is required for the interferometric measurement. 

LISA interferometry utilizes coherent optical transponders in order to compensate for the large diffraction losses caused 
by the immense distances between the spacecraft. A schematic representation of the LISA interferometer (restricted to 
two arms) is shown in Figure 1. Each of the three spacecraft, labeled contain two instrumental payloads 

labeled ‘A1 ‘A2’, ‘B1 etc. These payloads consist of laser, an optical bench, and a phase measurement system. The 
phase-measurement system measures the phase difference between the local laser and the laser on the adjacent bench as 
well asthe phase difference between the local laser and the received beam from the far spacecraft. These measurements 
are mani pul ated and combi ned to generate si gnal s that suppress I aser phase noi se whi I e retai ni ng the gravi tati onal wave 
signals. 

For example, consider the laser associated with payload A1 on spacecraft A as the ‘master’ laser. The back-link fiber 
measurement is used to measure the phase difference between lasers A2 and A1 . By adjusting the phase of laser A2 to 
track that of I aser A 1 (a phase-1 ocked I oop) , the two beams travel i ng to the di stant spacecraft behave as i f they ori gi nated 
from the same laser. Similarly, lasers B1 and Cl track the incoming beam from the spacecraft A. With all of these loops 
in operation, the two remaining phase measurements on spacecraft A produce the foil owing signals: 


<t>Aj(t ) = (p(t) - 4>(t - t aj ) ;j = B,C 


W here 0 (t) i s the phase of the master I aser and r AJ i s the round-tri p I i ght travel ti me between spacecraft A and j . 
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Each of these signals contains the desired strain information but is dominated by the random fluctuations of the laser 
phase. A technique known as Time Delay I nterferometry (TDI) [6,7] combines the two signals with varying time delays 
in order to suppress these fluctuations while retaining the gravitational wave signal. The noise suppression of TDI is 
primarily limited by knowledge of the nominal arm lengths, t A j -. With the assumption that a ranging system can track 
these arm lengths to 1 m, TDI can meet the LISA measurement requirements if the phase fluctuations of the master laser 
are below 


282 Hz 
VHz 



O.lmHz < / < 100 mHz 


(1) 


The expected free-running noise of the LISA lasers is 10kHz/Vtfz(^r), hence it is necessary to actively stabilize the 

laser frequency in order to meet the above requirement. Several methods of frequency stabilization have been proposed, 
including traditional laboratory techniques such as spectroscopic references [8] and high-finesse optical cavities [9]. 
Here we consider a technique known as arm-locking [10-17], which uses the arms of the LISA constellation as a 
frequency reference. 

The basic concept of arm-locking is to use the information contained in the phasemeter signals <p Aj to estimate the 
phase'frequency noise of the master laser and thereby correct it. This concept was first proposed by Sheard, et al. [10] 
and has subsequently been further refined and studied by others [11-17]. In this work, we apply the techniques of 
classical optimal control theory to this problem. 

ARM LOCKING 

For this initial study, the phase-locked loops for the local and remote lasers are not modeled. This assumption should not 
impact the developments presented herein as their contribution is well above the LISA measurement band. Additional 
noise sources that enter the phasemeter measurements on spacecraft A are considered, leading to a refinement of the 
laser phase dynamics. 

• v laser denotes the uncompensated frequency noise of the master laser. 

• v$f represents the frequency noise induced by the motion of the spacecraft A with respect to spacecraft B along 
theAB arm. 

• Vsc represents the frequency noise i nduced by the motion of the spacecraft A with respect to spacecraft C along 
the AC arm. 

• v ciock denotes the frequency noi se brought on due to the uncertai nty i n the cl ock i n the A B phasemeter readout. 

• v clock denotes the frequency noi se brought on due to the uncertai nty i n the cl ock i n the AC phasemeter readout. 

• ^shot represents the shot noise in theAB phasemeter readout. 

• v sh 0t represents the shot noise in the AC phasemeter readout. 


The parameter “ALS” refers to the arm locking sensor, which is designed to recover the optical phase within the LISA 
bandwidth using the outputs of the phasemeters The parameter “K(s)” denotes the compensator for the arm locking 
loop. Theinternal dynamics of the phasemeter, the actuation system for laser frequency control, and the phase-locked 
loops at the local and remote spacecraft, are not considered in the paper, as they do not have any significant impact to the 


2 



design or performance of the A LS and the compensator in the LISA band. Here, the dynamics of the optical signalsis 
expressed in terms of frequency for convenience. The outputs of the phasemeters are assumed to be processed to act as a 
frequency-meter. The outputs of the phasemeters are given as 


Pab(s) = (1 - e~ STAB ')v 1 (s) + vf ock (s) + V^ ot (s) + (s) (2) 


Pac(s) = (1 - e- ST ^) Vl (s) + V^ ock (s) + V^ ot (s) + v££(s) (3) 


TheALSoutput blends the signals from the phasemeters, i.e., 

S(s) = [SabCs) $«c(s)] {£"£>} 


Where S AB (s) and Sac 0) represent the elements of the ALS corresponding to the AB and AC phasemeters, respectively. 
The closed-loop system dynamics may now be written as 

[1 + K(s)G(s)] Vl (s) = -K (s)S AB (s)Vnoise (S) - K(s)S AC (s)V^ ise (S) + V Zaser (S) (5) 


Where 

G(s) = S AB (s)L ab (s) + S ac (s)L ac (s') 

Lab(s) = (1 - e~ STAB ) 

Lac(s) = (1 - e - "") 

VSStoCs) - v$* ck (s) + vZt (s) + Vsc(s) 

y noised) = v clocks) + v shot 0) + v£ C (s) 

The term [1+K(s)G(s)] isthe return difference, and K(s)G(s) denotes the loop gain for the arm locking loop. The laser 
noise with arm locking loop closed is written as 

ViO) = -[1 + K (s) G (s)] -1 ff (s) [S^ 0) V^ se (s) + 5' xc (s)v^j se (5)] (6) 

+ [1 + /f(s)G(s)]" 1 V Jajer (S) 


The challenge isto design an ALS, S(s), and acompensator, K(s), which would stabilize v t (s) to below required levels 
in the LISA band. 
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NOISE MODELS 


Laser Noise 

Arm locking can be implemented for afree running laser, or a laser that is pre-stabilized to a local reference, e.g., cavity. 
Although the latter option provides a more benign noise environment for arm locking, it adds additional complexities, 
such as hardware, mass, cost, etc. Moreover, si nee one of the objectives of the paper is to assess the feasibility of arm 
locking to meet LISA’s laser noise stability requirements, the laser is assumed to be free running, with no pre- 
stabilization. This represents a more challenging case for arm locking. The laser frequency noise model is based on 
laboratory measurements of afree running Nd:YAG NPRO laser. The spectrum of the noise drops with a slope of 1/f, 
and is provided below. 


V laser kHz 



Shot Noise 

Shot noise enters as additive noise at each phasemeter measurement. It is due to the quantization of detected photonsin 
the photoreciever and affects both the amplitude and phase of the signal. For an interference of two beams with powers 
P, and P 2 , assuming perfect contrast, the shot noise will affect the phase as 


0shot — 



WhereA is Planck’s constant, c isthe speed of light, tj isthe quantum efficiency of the detector, and A is the wavelength 
of the light. Putting in the constants for our 1 micron light and assuming rj= /thelevelis 


4> S hot — 69nCycles/'fHz 

Thisnoisewill be white in Fourier frequencies up to optical frequencies where it will roll-off. The statistics will be 
somewhat compl i cated but we can assume that we have enough photons that a Gaussian approxi mati on i s sufficient. T o 
convert this to afrequency noise in Hz, we simply multiply by 2*pi*f 



v S hot = 0.43 fiHz 




For LISA, we can assume 100 pW of light received from thefar SC is interfered with alocal beam of ~1 mW. Thisgives 
afrequency shot noise of 43^Hz/rtHz*(f/1Hz). 
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Clock Noise 


Clock noi9e enters in the phasemeter. This can affect both the PLLsat the far SC and the arm-locking sensor. A rough 
estimate of the magnitude of the clock noise is 


Vclock 





Where v beat is the frequency of the beat note. Here, a range of 5-1 5 M Hz is used for v beat . The clock noise will be 
independent on each spacecraft. For arm-locking it should be correlated between the two phasemeter signals at the corner 
spacecraft. 


Spacecraft Motion 

The current LISA baseline cal Is for using the bench-to-bench measurement to form the arm-locking error signals. This 
means that the spacecraft motion will be present in the sensors. This noise should be added as a varying time-del ay 
between the spacecraft. In other words, it affects the transfer function of the arm rather than the noise. However, given its 
small size (~2.5nm/rtHz) it is reasonable to treat it as an additive noise term in the phasemeter signal. The noise can be 
approximated by the Doppler effects of spacecraft motion as perceived in the phasemeter. 


With A = 1.064 jum, 


v sc 


-© 


2.5 X 10 ~ 9 Hz, 




v sc — 14.8 mHz 


■^Ge) 


ARM LOCKING SENSOR 

The arm locking sensor (ALS) isafilter or software which uses the measurement data from the AB and AC phasemeters 
to provide an estimate of the current laser frequency/phase, with some acceptable latency. This estimate is used by the 
dual arm locking compensator to correct for frequency/phase errors The approach taken here for the design of ALS is 
based on optimal estimation theory. Assuming a discrete-time implementation of the estimation scheme, let At denote 
the sampling time. Then, the phasemerter dynamics provided in Eqs. (2) and (3) may be approximated in the z-domain 
by 


Pab(z) = (1 -z fc ^8)v 1 (z) + V noise (^0 


Pac 0) = (1 - z~ k Ac)v 1 (z) + v£u. 00 (8) 
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Where k AB = ceil(r AB / At) and k AC = ceil(r AC / At). Without any loss in generality, assume k AB > k AC . Note that the 
discrete representations in Eqs. (7) and (8) should not impact laser frequency/phase estimation in the critical frequency 
band, which lies at or around the LISA band, if At is chosen small enough such that 1/ At is sufficiently outside the band 
of interest. With the phasemeter dynamics defined in the above, the dynamics of the arm locking loop may be written in 
state-space form 


x k + 1 =Ax k +Bu k + Tw k 


( 9 ) 


Vk = C x k + D u k + H w k + v k 


(10) 


Where the state and measurement vectors are defined as 


x k = pi* v i*-i - v i H- kAB y 

y k - { P AB k P AC k } T 

The state vector includes the laser frequency at the current cycle, the laser frequency from the previous cycle, and so on, 
all the way to the laser frequency from fc„ B cycles ago. However, note that the current cycle is delayed from the current 
time by one cycle, as it does not include the instantaneous impact of laser frequency actuation and laser noise. The state 
matri x A i s sparse with non-zero el ements be) ow the di agonal , i.e.. 
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Note the non-zero element in the first row of the output influence matrix, C, is the last element, corresponding to the 
state at /recycles ago, while the non-zero element in the second row is element no. k AC . Finally, the process noise 
vector, w k , and the measurement noise vector, v k , are defined as 


Wfc — V laser k 


v k = {v 


AB 


..Ac 


noise fc ''noisefc) 
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The ALS must provide an estimate of the laser frequency to be used in the arm locking controls loop. To accomplish 
this, optimal filtering is proposed. In particular, since only steady-state processes are involved and/or are of main 
concern, Steady-state Kalman filtering is used to compute an estimate of the laser frequency [18-19]. The steady-state 
Kalman filter may be realized in state-space form and is represented by the following equations. 


x k +i (-) = A x k (-) + Bu k + L (y k - Cx fc (-) - Dit fc ) 


( 11 ) 


*fc(+) = **(-) + M(y fc - Cx fe (-) - Du fe ) 


(12) 


9k - C £fc(+) + Du k 


(13) 


Here the vector 2 k (~) represents an optimal estimate of the state vector x k based on all prior measurements, while the 
vector x k (+) represents an optimal estimate which includes the current measurement as well. Equation(ll) provides the 
time update, while Eqs. (12) and (13) provide the measurement updates. The only estimate that isof concern for ALS is 
the most current estimate of the laser frequency the Kalman filter can provide. And that estimate is the first element of 
the vector 2 k (+). 


Filter Realization 

It can be shown that due to the sparse structure of the plant, represented by Eq. (9) and (10), the implementation of the 
optimal ALS (OALS) filter may be separated from the control input, such that the OALS filter assumes the following 
from. 


*k+i(-) = A %(-) + L (y k - Cx fe (-) ) 


(14) 


Vi fc (+) = v lfc (-) + Mi(y ft - Cx fc (- )) 


(15) 


In other words, one can show that for every controller and filter combination for the system in Eqs. (11 )-(1 3) there is an 
equivalent filter and controller combination, with a filter represented by Eqs. (14) and (15). Of course, one can continue 
to use the filter equations (11)-(13), which includes the control signal, or a combination of the two implementation 
options. 

The resulting OALS filter, represented by Eqs. (11)-(13), or Eqs. (14)-(15), can be high order, depending on the sample 
time At. However, the system matrices of the ALS filter are sparse, and hence the algebraic propagation can be done 
quite efficiently for flight software implementation using sparse matrix operations. Moreover, the dynamic range of 
interest resides primarily at or around the LISA measurement band. As such, a balanced reduction procedure may be 
followed to obtain a low-order approximation of the OALS filter for implementation and analysis. Reference [20] 
provides a thorough description of balanced realization and reduction process. 
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COMPENSATOR DESIGN 


The arm locking controller U9es the estimated frequency error provided by the OALS, save any commanded offset, to 
adjust thefrequency of the laser. The laser actuation dynamics, which may include thermal and piezo actuations, is not 
included in the compensator design or follow-up analysis. The objective of the compensator design isto suppress the 
laser frequency noise to alevel below pre-TDI requirement (Eq. (1)), whileensuring that the sensing-related noise, e.g., 
clock noise, shot noise, spacecraft motion noise, does not amplify unreasonably outside the bandwidth of interest. 
Furthermore, the arm locking compensator has to contend with and suppress spikes in the plant dynamics, within the 
L I SA band, that are brought on by the zeros due to round-tri p delays i n the two arms. The compensator has to provide at 
I east 6dB of gai n margi n and 30° of phase margi n for the arm I ocki ng I oop. The arm I ocki ng compensators were 
designed using the classical frequency domain approach. The designed compensators generally include a low-frequency 
lead-lag/lag component, followed by afirst-order attenuation filter. 

ANALYSIS RESULTS AND DISCUSSION 

The proposed arm locking control strategy was applied to a LISA constellation configuration. Here, the round-trip times 
for the laser signals were assumed to be nominally at 33 and 32.4 seconds, i.e., t ab = 33 and t ac = 32.4. The optimal 
arm locking sensor was designed using the optimal estimation approach described earlier. The sample time for filter was 
chosen at 0.1 sec, leading to afilter order of 331 . The noise models used for the ALS design were similar to those 
described earlier. The models were an exact match within the LISA band. However, below and above band the emphasis 
on some models were changed to provide desired loop shaping. Thefilter implementation was based on a parti al control 
signal feedback (see Eqs. (1 1)-(13)). A compensator was designed, whose Bode plot is shown in Fig. 2. As expected, the 
highest gain is concentrated below and within the LISA band. ThisOALSfilter and compensator combination is 
referred to as controller K1. Figure 3 shows the I oop gain, with the attained stability margins. Both gain and phase 
margins meet the requirements Laser frequency noise suppression is demonstrated in Fig. 4. It is observed that the laser 
frequency noise is well below 282Hz/-fHz throughout the LISA band, and hence meets LISA’s pre-TDI requirements 
The response is dominated mostly by the free-running laser frequency noise. The contributions from the clock noise and 
spacecraft jitter provideafloorfor noise suppression. Note that the sequence of peaks within and close to the LISA band 
occur at the frequencies equal to multi pies of the inverse round-trip delay of the arms These peaks have been effectively 
suppressed by using the estimated laser frequency from OALS. This is made possible by using phased requency 
information from two arms If one had to resort to single arm locking, then these peaks are not amenable to suppression 
as they lie on the zeros of the plant dynamics at multiples of the inverse round-trip delay. 


ALS Sensitivity 

The arm locking sensor filter isbased on discrete-time steps, which means that round-trip time for each arm, and 
therefore the difference between them, is captured at that quantization level . Since practical considerations, such as flight 
software implementation, limit the fineness of the step-size, it is imperative to assure the OALS is robust to errorsinthe 
representati on of the round-tri p ti mes, and hence the di f f erence between them. 

The closed-loop arm locking analysis was repeated with controller K1, which was designed based on t ab = 33 and 
r AC = 32.4. However, for the analysis-r^ was set to 32.5, resulting in the difference in round-trip times of 0.5 seconds. 
The performance of arm locking controller is provided in Fig. 5, where little or no change is observed in the level of laser 
frequency suppression from the nominal case. Hence, the design appears to be insensitive to representation errors in the 
round-trip time for a single arm, or the difference in the round-trip times between the two arms, of O.lsecor smaller. On- 
orbit implementation would switch the OALS filter and compensators for their most optimal at a time quantization of 
equal or better than 0.1 seconds. Therefore, the proposed OALSfilter and compensation approach is feasible for on-orbit 
implementation. 

In order to exercise the range of feasibility of the OALS filter, anew analysis was performed with controller K1. 
However, thistimer^c was set to 32.95, resulting in the difference in round-trip times of 0.05 seconds. Figure6 
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illustrates the laser noise frequency suppression for this setup. Although LISA pre-TDI requirements are met for most of 
the LISA band, however significant poke-throughs are observed at or near the multi pies of inverse round-trip times. 

Next, anew OALSfilter and compensator were designed using the best representation of the round-trip timesfor the 
arms. This representation set the values of r^and z xc to 33 sec. and 32.9 sec., respectively. The performance of the 
optimal filter isshown in Fig. 7. It is observed that the laser frequency noise is suppressed well below 282 Hz/fHz 
throughout the LISA band, and hence meets LISA’s pre-TDI requirements. 

As the differential arm length approaches zero, the ripple in the closed-loop OALS system near the 1/t frequencies 
increases and eventually becomes larger than the free-running noise. Thisisbecausethenullsintheindividual arm 
responses are degenerate and no combination of the two can be used to measure phase noise at these Fourier frequencies. 
This is a common feature of all loops that utilize two LISA arms [15]. In the nominal case when all three LISA arms are 
available, it will be possible to select the master SC to ensure that the differential arm-length is sufficiently great to avoid 
this problem. However, in the event of a loss of one of the arms, the master SC isfixed and there will be times when the 
differential arm-length goes to zero. 

By systematically varying the differential arm-length, it wasfound that the OALSfilter presented above will meet the 
pre-TDI frequency noise requirements so long as the differential arm-length exceeds 2msec. Thefraction of the orbit 
with differential arm-lengths smaller than this varies depending on the arm that is lost and the exact orbital model. Using 
asimplified model (first-order in eccentricity, no perturbations from outer planets, etc.), it isesti mated that a differential 
arm-length of I ess than 2msec will last for a period of ~1/2day. Therefore, inthecasethat alink is lost, theOALSfilter 
will still satisfy mission requirements for 99.9 % of the mission duration. 


CONCLUDING REMARKS 

The problem of arm-locking for LISA using phasemeter/frequency-meter datafrom two of its arms was considered. 
Using classical optimal estimation theory, an approach for designing optimal filters to recover the laser frequency noise 
within the LISA band was presented. Feasible compensators for the arm locking loop were designed and implemented. . 
Numerical analysis demonstrated that the proposed strategy to use optimal estimation of laser phased requency for arm 
locking applications is feasible. It is robust to variations in representation errorsfor round-trip times of the arms, and has 
the advantage of not requiring precise knowledge of the round-trip times or their difference. 

An important effect that has recently been recognized for arm-locking isthat of Doppler estimation error in the 
phasemeter. Variations in the relative velocity between the spacecraft add atime-varying Doppler shift to the 
phasemeter signal. This component is subtracted from the phasemeter signal using an orbital model. Thismodel isnot 
perfect however and some residual Doppler shifts remain and affect both the residual noise and dynamic range of the 
laser. This effect will be considered in subsequent studies. 
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Figure 1 Block Diagram for Arm Locking 


Soft Diagram 
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Figure 2 Bode Plot for the Arm Locking Compensator 
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Figure 3 Loop Gain for the Arm Locking Loop 



Figure 4 PSD for Laser Noise Frequency: Controller K1 
= 33 and = 32.4 




Figure 5 PSD for Laser Noise Frequency: Controller Kl, Figure 6 PSD for Laser Noise Frequency: Controller Kl, 

x A g — 33 and x A q — 32. 5 ^ab = 33 and = 32. 95 



Figure 7 PSD for Laser Noise Frequency: Optimal Controller, 
x AB = 33 and x AC = 32. 95 
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